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Abstract 

We study the evolution of the Weyl curvature invariant in all spatially 
homogeneous universe models containing a non-tilted 7- law perfect fluid. 
We investigate all the Bianchi and Thurston type universe models and 
calculate the asymptotic evolution of Weyl curvature invariant for generic 
solutions to the Einstein field equations. The influence of compact topol- 
ogy on Bianchi types with hyperbolic space sections is also considered. 
Special emphasis is placed on the late-time behaviour where several in- 
teresting properties of the Weyl curvature invariant occur. The late-time 
behaviour is classified into five distinctive categories. It is found that for a 
large class of models, the generic late-time behaviour the Weyl curvature 
invariant is to dominate the Ricci invariant at late times. This behaviour 
occurs in universe models which have future attractors that are plane- 
wave spacetimes, for which all scalar curvature invariants vanish. The 
overall behaviour of the Weyl curvature invariant is discussed in relation 
to the proposal that some function of the Weyl tensor or its invariants 
should play the role of a gravitational 'entropy' for cosmological evolu- 
tion. In particular, it is found that for all ever-expanding models the 
measure of gravitational entropy proposed by Gr0n and Hervik increases 
at late times. 



1 Introduction 

What is the entropy of a gravitational field? Following the insights of Hawking 
and Bekenstein jj], |^], an answer has been only been provided for a subset of 
situations characterised by static and stationary spacetimes with event hori- 
zons. In particular, black holes are found to be black bodies. They possess a 
well-defined entropy. They obey the laws of equilibrium thermodynamics and, 
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when perturbed, their fluctuations obey the laws of close-to-equilibrium thermo- 
dynamics ||. Yet, despite these developments it remains an unsolved problem 
whether an analogous gravitational entropy can be defined for non-stationary 
spacetimes of the sort that characterise most relativistic cosmologies that are 
considered as models for the past, present, and future of our visible universe. 
Unfortunately, a solution of this problem appears to involve the solution of a 
number of separate hard problems in gravitation physics. It requires a rigor- 
ous measure of the probability of different cosmological initial conditions to be 
found. It also requires us to be sure that we have accounted for all the possible 
contributions to the gravitational entropy: there may be quantum, classical, 
geometrical, topological, thermal, and dynamical contributions that need to be 
quantified and added 0] . The nature of the states of zero and maximum entropy 
also need to be defined and the influence of possible events like inflation must be 
fully accounted for entropically, and all must be suitably coordinate covariant. 
It has also been argued [|[ ^| that the whole concept of the arrow of time and 
thermodynamics relies upon the concept of a gravitational entropy, and hence 
leads all the way back to fundamental questions about the origin of time and 
the nature of the initial state of the universe. 

Gravitation has a tendency to amplify small inhomogeneities in the distribu- 
tion of matter. Yet, our present observations of the spectrum and temperature 
isotropy of the microwave background radiation show that the radiation was 
very close to thermal equilibrium when it was last-scattered, and that fluctua- 
tions from the Friedmann-Robertson- Walker (FRW) spacetime metric are only 
of order 10~ 5 . If there does exist a gravitational entropy that grows with time 
to reflect the action of gravitational instability then it should have the property 
that it vanishes (or at least is small) for a homogeneous and isotropic metric and 
be maximal for a black hole of the same total mass. This type of observation has 
led some investigators to propose that any measure of cosmological gravitational 
entropy should be proportional to the deviation of the universe from the FRW 
model [Q. Hence, the present gravitational entropy of the universe is in some 
sense small and so must have been even smaller in the past, perhaps providing 
some rationale for assuming the initial state of the universe to have zero entropy 
and hence be FRW up to minimal quantum gravitational fluctuations. 

Based on this observation, Penrose fl suggested that the Weyl curvature 
invariant, given by 

C a ^ s C afilS = R a ^ S R a07S - 2R al3 R al3 + |i? 2 , (1) 

should be very small near the initial singularity, and then grow thereafter behav- 
ing as a proper measure of the gravitational entropy. This in all its editions, is 
what we will call the Weyl Curvature Conjecture (WCC). Wainwright and An- 
derson H had a different version of the WCC. They expressed the conjecture 
in terms of the relative magnitude of the Weyl and Ricci curvature invariants: 

However, based on observations from two different universe models and taking 
into account that the entropy should scale with the volume, Gr0n and Hervik 
[ p)| , suggested that a better measure of a gravitational entropy would be 

s = VW\p- (3) 
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Most other studies of the WCC @ [l| |T§ |T|, [[?], [l| use the square of 
the Weyl scalar, C a ^ jS C a i3 1 s, or P, as the putative gravitational entropy for 
cosmological models. However, these candidates do not scale as the volume 
and so they do not capture correctly the dependence of a Weyl entropy on the 
volume. 

We introduce the expansion-normalised shear: 

_ 3 
h ~ 2 02 • 

The shear tensor is a symmetric and trace-free tensor defined by 

V»u = Qfiu - ^(g^ + u^u v ). 

Here, u M is the normalised four-velocity of the comoving fluid; = u^ v is the 
expansion tensor of the comoving fluid and 9 = 0^ the volume expansion rate. 
The shear tensor tells us how anisotropic the universe is; in particular, £ = if 
and only if the universe is isotropic. From the generalised Friedmann equation 
(8ttG = c = 1) 

\0* = \a>"v» v -\WR+p + K. (4) 

Since the density, /?, is positive, £ will be bounded by 

< £ < 1 (5) 

whenever the spatial Ricci-scalar is non-positive: < 0. This inequality will 
be fulfilled in all of the homogeneous models, except in the type IX and the 
Kantowski-Sachs closed universes. 

One of the aims of this paper will be to investigate the quantity S for all 
homogeneous models, and determine whether or not it has the appropriate be- 
haviour for a gravitational entropy. Hence, we will try to determine whether 
or not S increases during the evolution of all homogeneous universes. We will 
also investigate the quantity P, because this quantity tells us something of the 
physical properties of the universe model under consideration and the general 
properties of the evolution of anisotropic universes. We will discuss, unless 
stated otherwise, only spatially homogeneous cosmologies with a non-tilted 7- 
law perfect fluid. We also assume that the matter sources obey the strong energy 
condition; the effects from inflationary fluids are considered elsewhere jflj . Wc 
will investigate all the possible spatially homogeneous models. Of special inter- 
est are the eight Thurston geometries jll| which are related to the classification 
of topologies in three dimensions [f20[. T hese models correspond to the Bianchi 
types in a rather interesting way ]21|, [22|, |23] | . The possibility of compactification 
of these Thurston models has induced a renewed interest in the homogeneous 
models of Thurston and Bianchi type ^5|, |2a 2?J . 

This paper is organised as follows. In section |2 we investigate all the models 
that have the FRW models as a special case. The remaining Thurston models 
are investigated in section |[ and in section |J we investigate the remaining 
Bianchi models of non-Thurston type. We summarise our results in the final 
section. 
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2 Models with maximally symmetric Thurston 
geometries 



First, we will consider models which contain the FRW universes as special cases. 
The spatial sections in these cases are either E 3 , H 3 , or S 3 . In the isotropic case 
these correspond to the Bianchi types I and VIIo, V and VII^, and IX, which 
contain the k — 0, — 1, and +1 FRW universes respectively. 



2.1 The case E 3 

For the flat case, ^ R = 0, we have already exhaustively examined the Bianchi 
type I models in earlier papers [n). An important class of solutions is 
provided by the Kasner vacuum solutions: 

ds 2 = -dt 2 + t 2 ^dx 2 + t 2 ^dy 2 + t 2 ^dz 2 

The exponents obey YliPi = YliPi = 1- The solution space of these solutions 
is a circle, and are usually referred to as the Kasner circle. IiQ pi = 1, then 
the line element is Taub's form of flat spacetime. This form of flat spacetime 
will be referred to as T, and lies on the Kasner circle. Note that T is one of 
two Locally Rotationally Symmetric (LRS) spacetimes on the Kasner circle. All 
vacuum Kasner solutions have maximal expansion-normalised shear: 

E = 1. 



We will assume a perfect fluid with equation of state linking the pressure p 
and density p : 

P= (7- l )P- 

For all 7 < 2 it can be shown that the Kasner circle is a past attractor^. 
Calculating the Weyl tensor, we find that near the initial singularity 

where t is the cosmological time. The quantities P and S defined in eqs. (||) 
and (Q), can be shown to behave as (assuming a Kasner-like behaviour) 

Pi oc f- 2 (6) 
Si oc f'- 1 

to leading order near the initial singularity. We note that P will diverge near the 
initial singularity for all 7 < 2; S will approach zero for 7 > 1 and diverge for 
7 < 1. If a fluid with 7 > 1 is assumed to dominate near the initial singularity 
then S will increase as t — * 0. Radiation fluids will dominate the early history 
of our universe so long as particle interaction and scattering times are shorter 
than the expansion time. This may not be the case for asymptotically-free 
interactions between 10 14 — 10 19 GeV but is likely to be restored (subject to 

lr The metric has permutation symmetry pj| , so any of the other exponents could equally 
well be 1. 

2 This is because the shear invariant cr^ v cr^ v , can be viewed as a 7 = 2 fluid in the Fried- 
mann equation mh. Hence, if the model contains a perfect fluid with 7 < 2, the shear will 
dominate the initial singularity. This need not to be the case when 7 = 2. 
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boundary effects on the long-wavelength part of the momentum distribution 
due to finite horizon effects and the breakdown of statistical mechanics) on 
approach to the Planck scale. Also, it is most likely that any self-interacting 
scalar fields present will become dominated by their kinetic energies as t — > 0. 
This is equivalent to domination by a 7 = 2 fluid. Thus an effective 7 > 1 in the 
neighbourhood of the initial singularity, is quite probable. It might not be the 
case if the potential energy of a scalar field dominated or scaled in proportion 
to its kinetic energy as t — > 0. An example of the latter type is provided by a 
scalar field ip with an exponential potential V((p) — Vq exp[— X(f\. 
We see that for the point T on the Kasner circle, we have 

{C a ^ 5 C aPlS ) T CX i-4. t 2(2- 7 ) 

Pt oc constant 

S T oc t (7) 

near the initial singularity. 

When we include a matter source with anisotropic pressure created by pho- 
tons or asymptotically-free particles in an anisotropic universe (since they are 
expected to be collisionless significantly below the Planck temperature) or a 
cosmic magnetic field, the behaviour near the initial singularity turns out to be 
quite complex. The behaviour of a Bianchi type I model with a magnetic field 
p9| or Yang-Mills field |30| is chaotic. Similar behaviour are found in Bianchi 
type IX and VIII. In the types VIII and IX, curvature effects cause the universe 
to "bounce" between different Kasner states; the magnetic field in the type I 
case, tends to mimic this effect. The Bianchi type I with pure magnetic field was 
investigated using the dynamical systems approach in ]29| . When the universe 
is in a Kasner epoch, the Weyl curvature invariant behaves as in the Kasner 
universe. The question to consider is: what happens to the Weyl curvature 
invariants during the bounces? 

As we approach the initial singularity, the duration of the bounces will be- 
come negligible compared to the duration of the Kasner epochs in between the 
bounces. We can write the Weyl curvature invariant as |3l| 

C a ^ s C a p^ = H*W 2 

where W is the dimensionless Weyl curvature invariant, and H = 9/3 is the 
Hubble factor. The quantity W 2 can be written so that it is only a polynomial 
(of degree 4) in the state space variables. The state space for Bianchi type I 
with a magnetic field is compact, and hence, all the variables are bounded and 
W 2 must therefore also be bounded at all times. Similarly, 

where 1Z 2 is a positive definite polynomial of the bounded state-space variables. 
When the universe is in a Kasner epoch, 1Z will vary as^| 

Uj = 0(t 2( - 2 -^) + o(<t(i-2cos0)^ 

3 When an expression like this contains two terms, we should only consider the term which 
is the dominant one in each situation. Sometimes the real expression contains cross-terms 
which will lie in-between the two. These terms will be omitted in this paper, mainly for 
brevity. 
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where </> is an angular variable on the Kasner circle. The first mode £ 2 ( 2_ t) comes 
from the perfect fluid term, while fsC 1 - 2008 ^) comes from the energy density of 
the magnetic field. The behaviour of P during these periods will therefore vary 
according to the value of the Kasner parameter </>. The quantity S will behave 
as 

S 7 = O (t 7 " 1 ) + O (t* ( 1+4 cos 0) ) 

during the Kasner epochs. If 7 > 1 this quantity will be increasing during 
this period. If 7 < 1 there will be Kasner epochs where S will decrease, and 
others where it will increase. On "average" , assuming that all the Kasner epochs 
occur with equal probability, the quantity S will increase. During a bounce, the 
magnetic field will experience a sharp peak in its energy density. This makes 1Z 
experience a peak as well. Hence, during this bounce, the quantity S will fall for 
a brief period, until it rises to a new Kasner level again. Thus, during Kasner 
oscillations, S will both decrease and increase in time. However, on average, S 
will increase continuously during these oscillations 

The Bianchi type VIIo universe model, the other Bianchi model which is 
compatible with E 3 spatial sections, has further interesting properties. For 
example, ifl<7<4/3 the solutions of this model will at late times isotropisc 
in terms of the shear, but not isotropise with respect to the Weyl tensor [g2[ 
It can be shown that if 7 < 4/3 

Zviio — > (8) 

while for 7 > 4/3 the shear will oscillate as t — ► 00 j34|. The Weyl tensor has 
a quite complicated behaviour; if 7 > 1 then W will oscillate with larger and 
larger amplitudes. The oscillations will diverge at late times in the sense thaf] 

limsup|W 2 | = 00 (9) 

t >oo 

This happens because the universe at late times will oscillate around an isotropic 
state. As the universe expands, the amplitude of the oscillations will be smaller 
and smaller, and hence, will isotropise with respect to the shear. But the os- 
cillations will become more and more rapid while W increases beyond bound 
§1- 

The scalar W measures the absolute value of Weyl curvature, while the 
quantity P measure the ratio of the "matter curvature" to the Weyl curvature. 
Strictly speaking, the measure P will oscillate as the quantity W does, but we 
can avoid this by either adopting the definition of the Weyl parameter as in 
[ p2[ [j3[ , or we can take the absolute value of W 2 and thereafter average over a 
complete period. We will do the latter, and denote these entities with a bar, for 
example P. In the case of VIIo we then get 

Pvii — > 00 j as t — > 00 

So in this case, even the "P-version" of the WCC holds. Since the universe is 
ever expanding, S will diverge at late times as well, as 

Svii — > 00, as t — ► 00. 

4 Note that in physically realistic universes, which begin their cosmic evolution at the Planck 
time, there is time f or v ery few oscillations to occur in the universe's history. 

5 The authors in j32| , |33| define a "Weyl parameter" W differently from that defined here. 
This makes their parameter increase monotonically at late times, while ours will oscillate. 
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2.2 The hyperbolic spaces H 3 

We will now consider the H 3 cases. They have ( 3 'R < 0. The future asymp- 
totics for these models are more easily obtained because they approach vacuum 
solutions in the future. All of the type V and VII^o solutions will, assuming 
that the matter content obeys the strong energy condition, evolve asymptot- 
ically towards plane- wave spacetimes |35|. For plane waves, all the curvature 
scalars are zero even though the metrics are anisotropic. In general, these open 
universe solutions will not have a zero Weyl tensor (except for a set of zero 
measure which correspond to the isotropic vacuum Milne universes). However, 
the late-time behaviour of VII/j models is to approach particular exact solutions 
first found by Lukash |3(| which are plane waves J|5|, Q .[] 

This behaviour indicates that the evolution of the Weyl scalar and S in this 
case is more subtle as we approach the future asymptote. The Weyl scalars will 
tend to zero even though the anisotropy freezes in with E approaching a constant 
value. Clearly, the Weyl scalar cannot capture the anisotropic plane- wave modes 
of the gravitational field. This can be considered as an objection to the argument 
that it is some function of the Weyl scalar which is the proper measure for the 
entropy. We need a more discriminating measure for the gravitational entropy, 
if it exists. It might also happen that these plane- wave spacetimes are isentropic 
modes corresponding to adiabatic perturbations of the gravitational field but 
we would need to know whether they are extendible into hermitian modes of a 
quantum cosmology. 

In the V and Vllh cases containing a perfect fluid, part of the Kasner circle 
will be a past attractor for all equation of state parameters 7 < 2. Hence, the 
discussion for the type I model near the initial singularity also applies to these 
models. 

As we already have mentioned, the future attractor is a Lukash plane-wave 
solution with E ^ if 2/3 < 7 < 2 and the open FRW model with E = if 
7 < 2/3. Using the exact solutions in the type V model with perfect fluid (see 
[pi]), we can show that when 2/3 < 7 < 2 the late-time behaviour for these 
solutions has a Weyl curvature invariant evolution 

{c a ^ 5 c aP , s ) v = t- 4 • [o{r*) + o(t-^- 2 )] . (10) 

The first of these terms comes from pure shear term effects in the absence of 
matter. The second term is driven by the matter interactions. In the presence of 
a perfect fluid, the pure shear term will always be sub-dominant for 2/3 < 7 < 2. 
Under the assumption 2/3 < 7 < 2 at late times we have 

P v oc t-i^ 

S v oc t&. (11) 

Even though this is derived from the exact solution of type V, one can show 
that eqs.(|io|) and ( pi] ) will also hold in the generic case (see Appendix A); the 
expression for Weyl curvature invariant in the type V case is remarkable simple. 
Thus, in the type V case, P is always decreasing and S is always increasing. 
Note also that this increase is due to the anisotropic modes; we assumed the 
universe to be homogeneous. 

6 There also exists an h = 4 / 11 Lukash vacuum solution found by Lukash j36| and a 7 = 2 
stiff fluid solution by Barrow |B7[. 
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The case VII/j is more difficult. Investigations by Barrow and Sonoda ]34|, p5||. 
and later by Hewitt and Wainwright |38|, |3l[] indicate that if2/3<7<2 then 
the generic solution asymptotes a plane-wave solution at late times. The plane- 
wave solutions have £ = constant and zero Weyl curvature invariant. The 
approach to the plane-wave asymptote, is subtle and a careful investigation of 
the Weyl curvature invariant in a neighbourhood of the plane- wave solutions is 
needed. The investigation and calculation is somewhat tedious, but the result 
is rather interesting: 



VII h 



oc t 



O t 



2(1+S + ) 



o t 



4E,+(3t-2) 



(121 



(see Appendix A on the nature of these asymptotes). Here, S+ is a constant 
parameter which tells us which plane- wave solution is approached as t — > oo, 
and E + obeys the inequality 



1 



(3 7 -2) < £+ < 0. 



The shear is related to £+ for these plane- wave solutions via the simple relation 



^-"plane-wave — | ^-f- 1 ■ 

Note the appearance of the shear-to-volume expansion rate parameter in the 
power of time rather than as an additive perturbation. This is a signal of the 
stability of the isotopic solution being decided at second (or higher) order be- 
cause the Lukash plane- waves are exact solutions of both the Einstein equations 
and the linearised Einstein equations. Thus a linearisation of the Einstein equa- 
tions about the Lukash solution gives a zero eigenvalue |3{|. Which of the two 
terms that in eq. (|l^) is dominant depends on the parameters £+ and 7. The 
Ricci scalar squared will in general decay as 



(R 



11 Rnv)vil h °^ t 



t 



,4S + + (3-y-2) 



Hence, 



vii h 



O t 



3(S + + (t-1)) 
-2S + 



4S + +(3 T -2) 
+ ( t 2 < 1 - 2S +) 



and the volume expands as \/\g\ oc t 3 ^ 1 2E +), so S will increase as 



= O [t 1 - 3E + +C t 2 < 1 - 2S +) 



Note that both P and S — > 00 as t — > 00. Both the Weyl invariant and the 
Ricci invariants approach zero, but the Weyl invariant does so at a slower rate 
than the Ricci invariant. Ricci has R^R^ oc p 2 , while the Weyl tensor contains 
shear modes as well. 

The work of Collins and Hawking |^| is of relevance here. This was a study of 
the stability of isotropic ever-expanding universes against spatially homogeneous 
perturbations of Bianchi type VII^ and was not restricted to perfect fluids and 
comoving fluids. In the case of VII^ perturbations to open FRW universes, 
isotropy was shown to be unstable. However, this well known result requires 



careful interpretation. The definition of stability used was asymptotic stability 
in the sense of Lyapunov; that is, as t — ► oo we need to have E — ► 0. However, in 
general we have found that these VII/j perturbations of the open FRW universe 
approach the Lukash plane-wave spacctimes which have E = constant (E — ► oo 
is forbidden by p > 0). Thus although isotropy is not asymptotically stable it is 
stable in the sense that the deviations from isotropy (E = 0) are always bounded 
[0 Hi- Note that the strong theorem of Collins and Hawking, assuming only 
the strong energy condition (7 > 2/3) is a reflection of the property that 'the 
matter doesn't matter' as t — > 00 in the open universe case. Thus, the behaviour 
determining the asymptotic stability is that of the vacuum solution and so does 
not have any significant impact upon why the present anisotropy is so low (see 
[[4l| for a discussion of how this may be explained). By contrast, in the case of 
type VIIo perturbations around the flat FRW model, the matter content of the 
universe is very influential and E decays to zero as t — > 00 only if the dominant 
form of matter is a perfect fluid with zero pressure (7 = 1) |59). Again, this 
is not very relevant to the present state of the universe since in the absence of 
inflation the bulk of its expansion e-folds since the Planck epoch have been in 
a radiation-dominated state for which isotropy is not stable in the Collins and 
Hawking sense. 

The space H 3 is infinite, and in the above cases we have assumed that we 
are dealing with an infinite and unbounded space H 3 . However, if we relax the 
assumption of global homogeneity and isotropy, we can get other spacctimes 
whose spatial section at each instant of time is a compact hyperbolic mani- 
fold (CHM). Interestingly, our earlier discussion must be completely revised in 
the compact case. The compact hyperbolic spaces are heavily constrained by 
Mostow's rigidity theorem; no anisotropic deformations of CHMs are allowed. 
All the Bianchi models of type V and VII^o with compact spatial sections have 
to be (locally) isotropic p3j . Hence, their Weyl tensors identically vanish.]] 

A more general investigation of the behaviour of S for CHMs is dependent 
on a calculation of the inhomogeneous eigenmodes of the CHM. This has not 
been done by date, but some numerical studies have been made |f43f . 

2.3 The case S 3 

Turning our attention to the closed universe case, the Bianchi type IX model 
has the possibility of chaotic dynamics. A closed universe will recollapse to a 
second singularity after a finite amount of proper time has elapsed if the matter 
obeys the strong energy condition, the positive pressure criteria, the dominant 
energy condition, and a matter regularity condition J44[ There can be an 
asymmetry between behaviour at the initial and final singularities which reflects 
the relative likelihood of different initial conditions. A special initial state that 
is chosen close to isotropy will deviate increasingly from isotropy on approach 

7 Inhomogeneous deformations are still allowed, but the spherical symmetry of H 3 is broken 
at the global scale when we go to its compact quotients. Thus, the inhomogeneous Tolman- 
Bondi model, which we discussed in refs. Jlo| must also be constrained. The matter 
content of the Tolman-Bondi models is determined by the dust-density p(r,t). If p is allowed 
to vary only as a function of the radius r inside a spherical shell of radius r;„j , where r;„j is the 
coordinate injectivity radius of the CHM, then the compactification can be done. Outside the 
shell of radius ri n j , the density p must be constant as a function of r so that the homogeneity 
can be "restored" outside this shell. The same arguments as for the general Tolman-Bondi 
model can now be made, and we get a S oc t 3 increase inside ri n j . 
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to any final singularity and exhibit chaotic spacetime oscillations induced by 
the anisotropy of the three-curvature. However, a quasi-isotropic initial state 
can be made general by the assumption that the matter content at the initial 
singularity is dominated by the influence of a 7 = 2 fluid or a scalar field 



with zero potential 37 . In this case the difference between the quantity W 
and P is also manifest. The quantity P is regular everywhere for the type IX 
case, except at the singularities. However, at the turning-point of the volume 
expansion H — 0, and so W will diverge there. 

When it comes to the initial singularity, the Bianchi type IX has a very 
peculiar behaviour. Near the initial singularity it displays chaotic behaviour in 
vacuum and with 7 ^ 2 perfect fluids (5], ||, ||, |g, || ||. There is the 
issue of the coordinate covariance of measures of chaos to be dealt with here, but 
this problem can be overcome |53). The geometrical picture is similar to that 
of the magnetic Bianchi type I universe which we have already discussed. The 
difference is in the Ricci tensor and in the matter terms. The Weyl curvature 
invariant will have the same property as in the magnetic Bianchi type I case. 
If the universe has a perfect fluid with equation of state parameter 7 then we 
would expect the same behaviour as in the Kasner universes. Hence, 

Six oc t 7 " 1 

during the Kasner epochs. The bounces induce downward spikes in <S, but on 
average we will have a steady increase of S with time until dust domination or 
potential dominated scalar- field evolution (7 < 1 effectively). 



3 Other models of Thurston type 

3.1 The case SX(2jk) 

Cosmological models with spatial sections 5L(2,R) can be invariant under two 
different Bianchi types, namely III and VIII. If we were to assume a type III 
invariance with 5L(2,R) spatial sections, we would find that this is only com- 
patible if the type III were part of a larger symmetry group. Hence, type III 
can be considered as a special case of VIII. 

Let us therefore investigate models of Bianchi type VIII. When we com- 
pactify the spatial space, this model becomes the most generic of the (locally) 
homogeneous models. The number of free parameters can be arbitrarily large as 
the topology of the space increases in complexity |24|, |2^, ^6) . In addition to this, 
the type VIII case is interesting for its dynamical behaviour as well. Not only 
does vacuum and perfect fluid (7 7^ 2) type VIII display chaotic behaviour near 
the initial singularity, as in the IX case, it can also display oscillatory behaviour 
at late times as in the VIIo evolution (54|, ^5). So, in this case, two rather pe- 
culiar features are manifested together, and by understanding the Bianchi type 
VIII could possibly be one of the keys for understanding the behaviour of S and 
the Weyl tensor in more generic models. 

Also the perfect fluid type VIII model has a similar behaviour as the type 
VIIq pQ] • For 1 < 7 < 2 the generic solution will approach the type III form of 
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flat space, but it does so in a way such that[] 

limsup|>V 2 | = oo. 

t >oo 

The shear has the asymptotic value E = \ in this case. Thus for the type VIII 
evolution we have 

Pviii~Pix, S V iii~Six as t — >0 
and for 1 < 7 < 2 (actually for 4/5 < 7 < 2) 

Wvin, Pviii, Svui — > 00 as t — ► 00. 
For 2/3 < 7 < 4/5, both W and E are finite and non-zero at late times; the 



solutions asymptote to the Collins type VI_i solution (see section 4.1). We refer 
to refs Q and [E7j for details. 

3.2 The case Sol 

Solvegeometry has the Bianchi Lie algebra VIo- The solutions will have the 
Collins VIo solution as a future attractor when 2/3 < 7 < 2. This solution is 
given by 

ds 2 = _ dt 2 + t 2 dj ,2 + t S=* , e 2c Xdy 2 + e -2c^2) 



where c = ^ 2 ^-f 1 ^ an< ^ P ~ ^-y'^ ^ ' ' ■ ^ ne shear is given explicitly by 

^ (37 -2) 2 



16 

We can also calculate the Weyl tensor for this spacetime: 



(r°(hSr \ - 2 (2- 7 )(3 7 -2) 2 (4- 5 7 ) 



The Ricci squared scalar is 
and hence, 



(2- 7 ) 2 (37 2 - 67 + 4) 

Via — ^,4^4 



2 ( 3 7- 2) 2 (4- 5 7 ) 

VI ° (2- 7 )(3 7 2 -6 7 + 4)- 

Unfortunately, P 2 < (the magnetic part of the Weyl curvature dominates the 
electric part) for 7 > 4/5, which makes P imaginary. Note also that P = 
for 7 = 4/5. Ignoring the sign (taking the absolute value) we can make P 2 
non-negative for all 7. In this case, we can define S: 

Svi oc t~i . 

Hence, since this is a future attractor solution, the generic solution will behave 
similarly at late times. 

The generic solution starts as a special Kasner solution, and thus P and S 
behave as in the Bianchi type I case near the initial singularity. 

8 Strictly speaking, these results are not rigorously proven; according to Wainwright they 
rely on a combination of stability analysis and numerical simulations |56|. 
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3.3 The case Nil 



Nil-geometry is invariant under the Bianchi type II group. The future attractor 
for generic 7-law perfect fluid solutions with 2/3<7<2is the Collins-Stewart 
solution: 

ds 2 = -dt 2 + t^r i^dx + -^ z dyj + (dy 2 + dz 2 ) 



where c = t/(2 - 7)^7 - 2) and p = (6 - 7)/(4 7 2 i 2 ). 

The shear and the Weyl curvature invariant can be readily calculated 

(3 7 -2) 2 
E// - 64 
(C a ^ 5 C a ^ 5 ) H oc r\ (13) 

Unfortunately, the Weyl proportionality factor can be negative, positive or zero. 
Taking the absolute value of the Weyl curvature invariant we obtain 

Pn = constant 

Sn oc tv. (14) 

The generic solution near the initial singularity will be approximated by a Kas- 
ner solution. 

3.4 The case ffxK 

This model corresponds to the Bianchi type III model. In the perfect fluid case, 
the nature of the future attractor depends on the equation of state parameter 
7. If 7 < 1, the future attractor is one of the Collins perfect fluid solutions. If 
7 > 1 then the generic solutions are future asymptotic to the type III plane- wave 
solution. Note also that the type III solutions can be considered as a special 
case of the VI^ solutions for h = — 1. The VI^ universes will in general be of 



non-Thurston type and will be considered in section 4.1. 

When 7 > 1, the solutions asymptote to a plane- wave solution and so we 
expect that the Weyl tensor decays faster than t~ A . Not many exact solutions are 
known in this case, but the LRS dust and radiation cases are known explicitly. 
Actually, in the LRS case with 7 > 1 we do not need the exact solutions to study 
the late-time behaviour. We can show that (see Appendix A)ifl< 7 <3/2, 
at late times the Weyl tensor will decay as 

(C^C^^^s^t-'-t-^l 

The late-time asymptotics still have shear and this is given by S = 1/4. The 
Ricci tensor will decay similarly, and so 

Piii.lrs = constant 

Shirrs oc t 2 . (15) 

The radiation case can be derived explicitly because the exact solutions are 
known. For 3/2 < 7 < 2 the Weyl tensor decays as 

[C af3 ^ 5 C a p 1 s) IU LRS oc t 4 • t 2 , 
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P and S vary as 



PlIIXRS <X t 2 ^- 3 

Sui, L rs « t 2 -<-\ (16) 

Hence, both P and S — > oo as t — * oo for 3/2 < 7 < 2. 

In the LRS dust case, 7 = 1, the Weyl curvature invariant will approach 
zero at late times as 

(C a ^ s C aPlS ) III LRS cx r A ■ {\nt)- 2 . 

The Ricci tensor also decreases at the same rate, so at late times P and S vary 
as 

Pin, lrs = constant 

Siii.lrs oc t 2 ln<. (17) 

The non-LRS case is slightly more troublesome. The type III dynamics 
can be considered as a special case of VI^ with h = — 1. The generic case is 
therefore considered in the VI/, section below. The LRS type III case is, as we 
explain in detail in Appendix A, a special case, and the behaviour of the Weyl 
curvature invariant is different in the general case than for the LRS type III 
case. If we want to compactify the type III geometry, we have to take an LRS 
model. The compactification procedure requires a higher symmetry than the 
three-dimensional type III symmetry. Hence, compact spatial geometry implies 
LRS. 

In the case 7 < 1, the future attr acto r will be the Collins VI ^ perfect fluid 
solutions with h = — 1 (see section |4.l| ). For this solution, the Weyl tensor 
decays as 

while the shear is constant: 

S/jj = |(3 T - 2) 2 . 
The Ricci square decays at the same rate, so 

Pin = constant, 

S ln oc t$. (18) 

Near the initial singularity, the solution is approximately that of a Kasncr 
solution. The LRS type III will approach the Taub point T on the Kasner circle. 
The analysis for the LRS type III shows that near the initial singularity, we get 
to leading order 

{G a ^C aPlS ) niLRS oc t^-t 2 ^ 
Pill, lrs = constant 

Slll,LRS oc t. (19) 

Hence, the type III behaves like the Bianchi type I universe near the initial 
singularity. 
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3.5 The case S 2 x R 



This includes the closed Kantowski-Sachs case and there is no corresponding 
Bianchi model. However, all the Kantowski-Sachs solutions are asymptotic to 
a plane symmetric Kasner model (not T) as t — > 0. Hence, near the initial 
singularity, we have 

The Ricci square is determined by the square of the energy density, and hence, 
as t — > 

Pks « f- 2 

Sks « f- 1 . (20) 

The Kantowski-Sachs models will recollapse for many types of fluids 0, 
Recent observations do suggest however, that we live in a state of accelerated ex- 
pansion. This type of behaviour can be explained by the late-time influence of a 
quintessence scalar field. In these cases, the KS universe can be ever-expanding. 
The KS models with a cosmological constant (7 = 0) contain an unusual unsta- 
ble particular solution which does not approach the de Sitter solution as t — > 00, 
[|58[ . The cosmic no hair theorem J59| is not violated because the 3-curvature is 
positive. This unusual behaviour has been used to argue for forms of dimension- 
ally selective inflation |6(J, however, the behaviour is an artefact of the S 2 x R 
symmetry and is unstable even within the KS class of solutions plj . 



4 Models of non-Thurston type 
4.1 The Bianchi type VI^. 

For this model, there exist future-attractor solutions for all 7 (as in the VIq 
case) . The nature of these attractors depends on the group parameter h (which 
is non-positive) and on the equation of state parameter 7. If 7 satisfies 

then the future attractor solution is the Collins VI^ perfect fluid solution with 
metric 

ds 2 = -dt 2 +t 2 dx 2 +t^Te rx2 ^r (t r ?e 2cx dy 2 +r r ?e- 2cx dz 2 ) , 



where c = >/(2 - 7) (37 - 2) and < r < 1 is related to h by h = -(2-^) 2 r 2 /c 2 . 
These solutions are universes with shear and some of their physical properties 
have been studied by one of us in connection with primordial nucleosynthesis 
[p2[ . The shear is found to be 

5Vz„ = ^(1-3/0(37 -2) 2 . 

The Weyl tensor in this case has a similar behaviour as in the type VIo case. 
For some values of 7 it is negative, and others positive. Ignoring the sign, as we 
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did in the VIq case, we get 



(C a ^ 5 C aPj5 ) VIh oc t- A 

Pvi h — constant 



Svi h 



t~ . 



If the equation of state parameter 7 fails to obey eq. (|T]) , and so 



h > 



7 



3 7 - 2' 



(22) 



(23) 



then the solutions appear to be asymptotic to plane-wave spacetimesQ. The 
plane-wave solutions have zero Weyl-curvature invariants, therefore we need to 
know how these solutions approach the future attractor. Unfortunately, only a 
few exact VI/j solutions with this behaviour are known. 

Nevertheless, we can show that when 7 obeys the inequality (|2^) , the be- 
haviour of a generic solution is 



{C a ^ s C aPlS ) VI <xt 



O t 



2(1 + B + ) 



O t 



4E + + (3-y-2) 



where E+ characterises the different plane- wave solutions (see Appendix A). 
The shear is directly related to this parameter via £ = The Ricci square 

goes as p 2 , and hence 



Pvi h =0[t 



3(S + +( T -1)) 



4S + + (3-,-2) 

+ 0[t *J=™+j 



Interestingly, in general P will grow unbounded as we approach the future at- 
tractor. Thus, in these solutions, the Ricci tensor will decrease faster than the 
shear modes captured by the Weyl curvature invariant. The Weyl tensor will 
decrease towards a plane-wave state, but as it does so, the Ricci tensor will be 
sub-dominant compared to the Weyl tensor. 

The volume expands as \f\g~\ oc t 3 /( 1_2S +) j thus iS will increase as 



& 



VIk 



3(S + +7) 

o ( t 1 - 2E + 



4S + +3 T 



(24) 



As mentioned earlier, these solutions are future attractor solutions, and 
hence, the generic solutions will have a similar behaviour for their Weyl ten- 
sors. Near the initial singularity, the generic behaviour is Kasner-like and the 
Weyl tensor behaves as in the Bianchi type I. 



4.2 The exceptional case 

The exceptional case of Bianchi type VI^jy g , has to be treated separately. The 
exceptional case has as many free parameters as the whole family VI/j with 
arbitrary h because the choice h — —1/9 results in the disappearance of some of 
the Einstein constraint equations. Unfortunately, the exceptional case has not 

9 This is not rigorously proven, but Hewitt and Wainwright pa, 3l| have conjectured that, 
except for sets of measure zero, these plane- wave solutions are global future attractors. The 
same is also the case for types IV and VII^. 
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been much studied in the literature. Stability properties were studied by Barrow 
and Sonoda [ p5| , a study of the consequences for primordial nucleosynthesis 
was made in J62| . and Hewitt J63[ , studied the evolution in a subspace of the 
whole dynamical system of VI*_]y g (namely the subspace with n a a = 0). In 
this subspace Hewitt and Barrow and Sonoda showed that for 2/3 < 7 < 10/9, 
7 = 10/9 and 10/9 < 7 < 2 the future attractors were the Collins Vl_ 1 / 9 perfect 
fluid solutions, Wainwright 7 = 10/9 solutions |6^] with metric 

ds 2 = -dt 2 + t 2 dx 2 + ti (e Ax dy + wtUx) 2 + th' 4Ax dz 2 

where w 2 = 9r 2 /4 — 1 and A — ry/6/5, and the Robinson- Trautman vacuum 
solution]^] with metric |3jJ 



ds 2 = -dt 2 + 7 4t 2 dx 2 + ^- e - x tUxdy + ^-e- 2x tidy 2 + e 4x tUz 2 , 
8 8 16 

respectively. We can readily calculate the Weyl tensors for these spacetimes. 
For 2/3 < 7 < 10/9 we get 

{C a ^ s C aMS ) VI , ocr\ 

v -1/9 

The Ricci tensor decays at the same rate, and hence 

Pvr oc constant 

v -1/9 

The Robinson- Trautman solution corresponds to a universe with shear £ = 
1/3, but has a zero Weyl curvature scalar. However, when 10/9 < 7 < 2 we can 
show that for any generic solution that approaches this solution, we have 

(C^ s C afflS ) VIli/g oct- 4 - [o(t"f) +o(t-^) 

(see Appendix B). Note that the two terms are equally dominant when 7 = 4/3 
(radiation). At late times the Ricci square decays faster towards zero than the 
Weyl curvature invariant: 

/ St— 11 \ / 9-1-10 \ 

P VIll/9 oc 0(t-r-)+0(rV-) 

S VIli/9 cx o{t 3 (^&) +o(t^ + &)) . (26) 

In this case the Weyl curvature invariant will dominate at late times. At early 
times, the generic solution approaches a certain point on the Kasner circle (not 
T), and hence, the analysis of Bianchi type I applies again. 

As mentioned, this analysis is just for a special subspace in the dynamical 
state space of VI* 1 , g . Nobody, at least to the authors' knowledge, has made 
a complete analysis of the whole dynamical behaviour of the exceptional case. 
However, page 171 in |H| appears to indicate that the exceptional case has a 
chaotic initial singularity. The future attractor seems to be a self-similar space- 
time, and hence the above results may hold in that case. Further investigation 
of this case is required. 

10 This solution is also sometimes called the Collinson- French solution. It contains no free 
parameters. 
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4.3 The Bianchi type IV 



It is conjectured that generically the Bianchi type IV solutions move asymptot- 
ically towards plane-wave solutions. The behaviour is similar to the VI^ case, 
and the late-time behaviour for the Weyl tensor is similar to the plane-wave 
case. Hence, we will have 



(C a ^ 5 C afilS ) IV oc i~ 4 • 



o 1 1 i< - +2B + ) ) + O I t 1 -•- 



(see Appendix A). Again the shear is given by S = |£+|. The Ricci square goes 
as p 2 , and hence 



VI, 



= o t 



3(S + + (t-1)) 



4S + + (3t-2) 

O I t 2 < 1 - 2S +) 



Finally, we have 



S 



vi h 



3(S + +t) 

O ( t 1 " 2S + ) +o I t 2(1 - 2S +» 



So, in this case the P and S — > oo at late times. Near the initial singularity, the 
type IV behaves like a Kasner universe, and hence, the Bianchi type I analysis 
applies again. 



5 Conclusion 

In this paper we have investigated the general behaviour of the Weyl curvature 
invariant for spatially homogeneous universes containing a perfect fluid as t — > 
and t — > oo. For all the models, except for the chaotic ones, the solutions 
asymptote in the past to points on the Kasner circle. The Kasner universe, 
which is a special vacuum case of the Bianchi type I model, has been treated 
exhaustively elsewhere [^7], |l(j. 

In the chaotic models, we saw that the Weyl curvature invariant C af3 ~ )S C a /3 1 s, 

the invariant P = — W i an d the invariant S = \J \g\P experience peaks 
and sinks. On the average, however, S increases steadily as the universe expands 
from the initial singularity. 

The late-time behaviour was found to fall into one of five categories: 

1. Ricci dominance. This category includes models for which P is decreas- 
ing at late times, indicating that the Ricci square dominates over the Weyl 
curvature invariant asymptotically. 

2. Weyl-Ricci balanced. In this category P is approximately constant at 
late times. The Ricci and Weyl tensors are approximately proportional at 
late times. 

3. Weyl dominance. In this category P — ► oo as t — * oo, but the expansion- 
normalized Weyl tensor, W 2 = C a,3lS C a p 1 sl H 4 , will decrease. Except in 
the exceptional VI^jy g case, these models evolve towards vacuum plane- 
wave spacetimes, but the Weyl curvature invariant decreases more slowly 
than the Ricci square. 
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Type 


Matter 


Category 


I 


2/3 < 7 < 2 


1 




2/3 < 7 < 4/3 + 7r^ 


1 




4/3 < 7 < 2+ mag.ficld 


2 


II 


2/3 < 7 < 2 


2 


HIlrs 


2/3 < 7 < 3/2 


2 




3/2 < 7 < 2 


3 


IV 


2/3 < 7 < 2 


3 


V 


2/3 < 7 < 2 


1 


VIo 


2/3 < 7 < 2 


2 


VL h 


on — n 

2/3<7<^ 


2 




^ < 7 < 2 


3 


VTLi/9 


2/3 < 7 < 10/9 


2 


10/9 < 7 < 2 


3 


VIIo 


2/3 < 7 < 1 


1 




1 < 7 < 2 


4 


VII,, 


2/3 < 7 < 2 


3 


VIII 


2/3 < 7 < 4/5 


2 




4/5 < 7 < 2 


4 


IX 


2/3 < 7 < 2 


5 


KS 


2/3 < 7 < 2 


5 


I-IX 


< 7 < 2/3 


1 



Table 1: Classification of the late-time behaviour. The matter content of the 
universe is labelled by the perfect fluid equation of state parameter, 7; 7r M „ 
indicates the presence of trace-free anisotropic stresses and type I also includes 
the case of a pure magnetic field. 



4. Extreme Weyl dominance. In these models, W — > 00 as t — * 00. 

5. Recollapsing models. These models have no late-time behaviour, they 
recollapse after a finite time subject to the usual technical energy condi- 
tions required for recollapse [Q. 

In the Table we have summarised our results, and in the rightmost column 
we have indicated to which of the five categories each of the spatially homo- 
geneous models belong. We note that category-3 cosmologies are not often 
discussed in the literature^ 1 ]. All of the plane-wave solutions of class B belong 
to this category. Note also that when we compactify the topology of the spatial 
sections , most of the models in this category disappear in the compact- 

ification procedure; types IV and Vlh cannot be compactified, type III must 
be locally rotationally symmetric (LRS), and VII^ must be isotropic. However, 
some of the type III spacetimes (namely those LRS with 7 > 3/2) in category 
3 can be compactified. 

A general criterion to decide when a dynamical system will end up in cate- 
gory 3 can be given. If the late-time asymptote is a vacuum spacetime which is 
not conformally flat, but nevertheless has a zero Weyl curvature invariant, then 

11 This behaviour seems to have been noted for a special class of the VI^ solutions in [|18[. 
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the quantity P will generically increase without bound. If, in addition, the state 
space of the corresponding dynamical system is also compact, then W has to be 
bounded (and thus decreasing), and hence, the dynamical system will end up in 
category 3. The reason for this is as follows. For a conformally flat spacetime, 
both the Weyl curvature invariant and the Weyl tensor itself vanish. Hence, 
any linear perturbation of this conformally flat solution will produce linear cor- 
rections to the components of the Weyl tensor. Since the linear terms are the 
lowest-order terms in the Weyl tensor, quadratic terms will be the lowest-order 
terms in the Weyl curvature invariant. For non-conformally flat spacetimes, the 
Weyl curvature tensor will not be zero. Hence, the linear terms will not be the 
lowest-order terms in the perturbation. Thus, in the Weyl curvature invariant 
some of the linear terms will survive even though the lowest-order (constant) 
terms could vanish identically (which happens for the plane- wave solutions). 
The electric and the magnetic parts of the Weyl tensor are not zero for plane- 
wave solutions, they just happen to be of equal magnitude and cancel when the 
invariant is formed. 

If the late-time asymptote is conformally flat, but not Ricci flat (hence, non- 
vacuum), then the system will typically end up in category 1. There is also 
a very small chance that the spacetime will end up in category 4 if the state 
space is non-compact. If the asymptote is both Ricci flat and conformally flat, 
all five categories are possible. If the state space is compact, then category 4 is 
ruled out. The fifth category depends to a certain extent on the geometry of the 
late-time attractor. For instance, inflationary spacetimes with p+3p < 0, are all 
governed by the "no hair" theorem which states that the late-time asymptote is 
a conformally flat non-vacuum spacetime. Hence, all ever-inflating cosmologies 
should end up in category 1. 

In all of the models investigated, except for those in category 5, S increases 
at late times. Hence, S behaves as expected by the WCC. It is not clear whether 
we should always expect a Weyl 'entropy' to increase if the initial conditions 
are close to some state of effective gravitational equilibrium. Of course, our 
study has only discussed spatially homogeneous universes and it remains to be 
seen whether small density and gravitational wave inhomogeneities evolve in 
a manner consistent with a thermodynamic interpretation of the Weyl curva- 
ture and whether the inclusion of significant inhomogeneities, black holes, and 
damping processes introduces significant new ingredients to the search for the 
elusive gravitational 'entropy'. 



Appendix A: The Weyl- Curvature Invariant for 
Bianchi types of class B. 

We will briefly show how the Weyl curvature invariant behaves near the plane- 
wave equilibrium points of type III, V, VI^, IV and VII/,,. The discussion is 
based on Hewitt and Wainwright's approach [[HJ |3^] (see also J|5| and 1 65 , |6^|). 



For all the plane-wave solutions we will have C a f3 1 sC- a ^' yS = 0, so the ques- 
tion we would like to answer is how the Weyl curvature invariant approaches 
zero at late times. The Weyl curvature invariant is given explicitly by 

C aPlS C a ^ s = 8ff 4 [{£ AB 8 AB + Q£\) - {H AB H AB + 6H 2 + )] 
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where the expansion-normalised electric parts and magnetic parts of the Weyl 
curvature are given by 



n+ 

c cAB 
CABO 



'Hab'H 



AB 



E+(1 + S H 
3A 



2 9 



-hA 

3 



6 



E(2E + - I) 2 +4l AT.+ + -(Nf- hA){A + (Nf. - hA) - 2E 



+4A 2 - 47V+A(2E+ - 1) + -hA{Nl - hA) 

o 

6[3Y, 2 + (Nl - hA) + 2%(Nl - hA) + tA + 6 A 2 
+12S+7V+A + AhAY] . 



(27) 



The shear is given by 



All the variables and their interpretations are explained in detail in refs. pl[ pq | . 
In particular, h is related to the Lie algebra parameter h via h = 1/h. The 
Bianchi types are related to h as follows. 



h< 
h = 
h>0 



IV 

vn h . 



(28) 



Note that type III is VI_i and type V is contained in the discussion as the 
special case of ft. = 0. 

In these variables the plane-wave equilibrium points are given by 



S = -S + (1 + S+) J A = 0, A=(l 
N\ = (l + E + )[A(l + E+)-3E+], 

where E + is a parameter which satisfies 

-K E+ < 
-1< E+ < 



s+) 2 , 



(29) 



h _ 
3-h ' 



if ft < 
if ft > 



(30) 



The eigenvalues of the dynamical system for these equilibrium points are given 

by 



{0, -4E+-(3 7 -2), -2[(1 + E+)±»JV + ]}, 



(31) 



where N+ is given by (29). The equilibrium points are future attractors when- 
ever 

1 

2). 



E + >--(3 7 - 

The second of these eigenvalues is directly related to the time evolution of 
the perfect fluid. If f2 is the expansion-normalised energy density, f2 will in a 
close neighbourhood of the plane-wave equilibrium points obey 



n' 



-4E, 



(3 7 -2)]0 
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-1 



Figure 1: In region A the eigenvalue — 4E + — (37 — 2) is larger than — 2(S + + 1). 
In region B, — 2(S + + 1) is the largest. 



where prime denotes derivation with respect to the time variable r defined by 

dt _ 1 
~aW ~ H' 

For the plane-wave solutions, we have O = 0. The dependence of the eigenvalues 
on the equation of state parameter is shown in Figure 1. 

Upon a small variation of the variables from these equilibrium points, the 
variation will decay as linear combination of factors exp(Air) where A^ are the 
eigenvalues (|3l|). Hence, we would naively expect the Weyl curvature to behave 
as 

W 2 = 0(exp(A l r)). 

But here we have to be careful. It can be shown that the generic behaviour of 
the Weyl curvature invariant is like this, but in some important cases it is not. 
For instance, for the radiative solution of LRS type III which we mentioned 



in section p.4| , all the linear terms from the variation around the plane-wave 
solution vanish. If we constrain ourself to lie in the invariant subspace A = 
£+iV + , 3X^ + KE = 0, A>0, then we would get 



W 2 = 



(1 + h) I (12ft + 7)Q + (69ft - 121)a+} + \Jl + ~h{27h + 5)v 

(32) 



(3 -ft) 

+ (quadratic terms in fl, cr + , v 



Here, f2 is a small variation of the density of the perfect fluid, and er+ and v 
are small variations of S + and N + respectively. The LRS type III is precisely 
this invariant subspace when h = — 1. Note that if and only if h = h = — 1, the 
linear terms in cr + and v will vanish identically. In the case were 7 < 3/2, f2 
will dominate at late times over ct+ and thus 

w 2 oc n 2 

at late times. Note also that R^R^ oc H 4 ■ il 2 which means that Pm ~ 
constant at late times for these LRS type III solutions. 
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The type V model is of special interest. The invariant subspace is given 
by h = 0, A > while E + = A = N + = 0. In this case the Weyl curvature 
invariant is particularly simple: 

Wy = 48(2E 2 + Eft). 

For completeness, we should mention that for all of the plane- wave solutions 
we have a simple inversion time evolution for the mean Hubble parameter. The 
simple Milne evolution arises only when the shear parameter S + = : 

H = r 1 . (33) 

1 - 2E+ ^ ; 

Appendix B: The Weyl- Curvature Invariant in 
the exceptional case of Bianchi type VI% 

9 

For the dynamical system of Hewitt ||| , the Weyl tensor is 

W 2 = -E 2 . (3 + 18S 2 2 + 4E 2 . - 167V 2 ) 

+32/V 2 (£+ + TV 2 - 3E 2 2 ) + 18S 2 2 (l + 3S 2 2 ) (34) 

E = ^ + 3E? a 
The Robinson -Trautman vacuum solution has 

E + ---, S ia _— , A^x-— ■ 

The eigenvalues of the linear system with respect to this equilibrium point are 
given by 

2 / /—A 10-9 7 



while the shear is 



3 V J 1 3 

In this case, the Weyl tensor will have linear terms in all of the three variables. 
Thus the Weyl tensor will decay as exp(A moa; r) where X max is the eigenvalue 
with the largest real part. The time variable t is related with the cosmological 
time via 

t oc e3 T 

near the equilibrium point. 
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